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Abstract 

Background: Aneuploidy, a karyotype deviating from multiples of a haploid chromosome set, affects the 
physiology of eukaryotes. In humans, aneuploidy is linked to pathological defects such as developmental 
abnormalities, mental retardation or cancer, but the underlying mechanisms remain elusive. There are many 
different types and origins of aneuploidy, but whether there is a uniform cellular response to aneuploidy in 
human cells has not been addressed so far. 

Results: Here we evaluate the transcription profiles of eleven trisomic and tetrasomic cell lines and two cell 
lines with complex aneuploid karyotypes. We identify a characteristic aneuploidy response pattern defined by 
upregulation of genes linked to endoplasmic reticulum, Golgi apparatus and lysosomes, and downregulation of 
DNA replication, transcription as well as ribosomes. Strikingly, complex aneuploidy elicits the same transcriptional 
changes as trisomy. To uncover the triggers of the response, we compared the profiles with transcription changes 
in human cells subjected to stress conditions. Interestingly, we found an overlap only with the response to 
treatment with the autophagy inhibitor bafilomycin Al. Finally, we identified 23 genes whose expression is 
significantly altered in all aneuploids and which may thus serve as aneuploidy markers. 

Conclusions: Our analysis shows that despite the variability in chromosome content, aneuploidy triggers uniform 
transcriptional response in human cells. A common response independent of the type of aneuploidy might be 
exploited as a novel target for cancer therapy. Moreover, the potential aneuploidy markers identified in our analysis 
might represent novel biomarkers to assess the malignant potential of a tumor. 



Background 

Aneuploidy, or a change in cellular chromosome numbers, 
has profound effects on the physiology of all eukaryotic 
cells analyzed to date [1]. Aneuploid yeasts are character- 
ized by slow growth, altered sensitivity to various stresses 
and increased genomic instability [2-4]. At the same time, 
aneuploidy drives genetic variability and cellular adap- 
tation capacity in yeast [5,6]. Plants are in general more 
tolerant to gene dosage changes, yet aneuploidy often 
impairs their vigour and alters their phenotype [7]. An- 
euploid mammals are rarely viable and the sporadic 
survivors are affected on multiple levels. In humans, 
aneuploidy is responsible for a substantial proportion 
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of spontaneous abortions and the rare survivors with 
trisomy of chromosome 13, 18 and 21 (Patau, Edward and 
Down syndrome, respectively) are severely handicapped; 
only trisomy 21 is compatible with survival until adult- 
hood [8]. Aneuploidy is also linked to cancer, as nearly 
90% of solid tumors and 75% of hematopoietic cancers 
show abnormal chromosome dosage [9]. Recently it has 
been shown that the occurrence of aneuploid cells in- 
creases with aging [10] and an increased incidence of 
aneuploidy in the brain has been linked to neurodegen- 
erative diseases [11]. 

The exact mechanisms underlying the detrimental 
effects of aneuploidy remain unclear, but it has been 
convincingly shown that they are caused by the expression 
of the extra genes on the supernumerary chromosome [3]. 
In most aneuploid cells the chromosome dosage changes 
lead to correlating changes in mRNA (e.g. [3,12-15]) as 
well as on protein levels [6,15,16]. These analyses further 
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revealed that besides the increased abundance of transcripts 
and proteins originating from the aneuploid chromosome, 
the expression of multiple other genes is altered as well 
[3,13,15,17]. This is likely a consequence of two different 
phenomena. First, an increase in gene copy number of a 
transcription factor or other regulatory factor might affect 
transcription levels of genes on other chromosomes [18]. 
Second, specific pathways might be activated in a cellular 
response to aneuploidy. Recent attempts to uncover the 
consequences of aneuploidy suggest that aneuploidy indeed 
instigates a specific response in eukaryotic cells [3,15,19]. 
Haploid yeast strains carrying an additional chromosome 
(hereafter referred to as disomes) exhibit a common tran- 
scriptional signature that has been previously identified in 
budding yeast as a so called environmental stress response 
that is triggered upon various exogenous stresses, such as 
oxidative stress, heat shock or slow growth [3,20]. A similar 
response was identified in a study comparing transcriptome 
data from disomic and complex aneuploid strains of 
budding yeast, partial aneuploids of fission yeast, aneu- 
ploid Arabidopsis thaliana plant cells, mouse cell lines 
with Robertsonian translocations that lead to trisomies, 
and human cells from patients with trisomy syndromes 
[19]. Additionally, comparative transcriptomics and prote- 
omics of model human trisomic and tetrasomic cells 
identified a common pattern in the transcriptional response 
to aneuploidy [15]. These studies pointed out similarities in 
the response to aneuploidy in most eukaryotes. At the same 
time, the results showed that the response in mammalian 
cells diverges from the response in other model organisms, 
because the correlation of the transcription changes in 
aneuploid mammalian cells with transcriptional changes 
in other aneuploid species is rather modest [19]. More- 
over, as there is no equivalent of the environmental stress 
response identified so far in mammalian cells, it remains 
unclear what triggers the transcriptional changes. 

To obtain a comprehensive insight into the changes 
specific for human cells, we evaluated in detail multiple 
tri- and tetrasomic cell lines as well as model complex 
aneuploid human cell lines with hypotetraploid karyo- 
types more similar to cancer cells. Remarkably, we found 
that the cellular responses to complex multichromoso- 
mal aneuploidy and trisomy or tetrasomy of one single 
chromosome closely correlate. As the identified pathway 
pattern resembles the cellular stress response, we com- 
pared the aneuploidy response to the transcriptional 
changes in human cells subjected to various stress con- 
ditions. Additionally, we have also identified 18 genes 
that were upregulated and 5 genes downregulated in all 
analyzed aneuploid cell lines and might thus serve as 
markers of aneuploidy. This is the first study which 
compares a variety of different human cell lines with 
aneuploidy of different types and origins. By uncovering 
a uniform aneuploidy response pattern our results outline 



the cellular consequences of an abnormal karyotype in 
human cells. 

Results 

Common cellular response to trisomy and tetrasomy in 
human cell lines 

Recently, we established de novo human aneuploid cell 
lines that were derived from diploid and chromosomally 
stable cells by a micronuclei-mediated chromosome trans- 
fer [15]. We performed transcriptional analysis of the 
original diploid cell lines and their trisomic and tetrasomic 
derivatives (HCT116: trisomy of chromosome 3 - 3/3; 
tetrasomy of chromosome 5 - 5/4; RPEl: trisomy of chro- 
mosome 5 and 12 - 5/3 12/3, trisomy of chromosome 
21 - 21/3, see Methods for further details) and calculated 
the aneuploid-to-diploid fold change for all detected 
mRNAs. Similarly, we calculated the relative aneuploid- 
to-diploid ratio to determine the transcriptional changes 
in trisomic human cells lines generated from the diploid 
colorectal adenocarcinoma cell line DLDl by introduction 
of an extra copy of chromosome 3, 7 or 13 [13]. To 
identify the similarities and differences among the 
transcriptional profiles of the model aneuploid cell 
lines in response to aneuploidy, we used an algorithm 
called 2-dimensional annotation analysis that quantifies 
the relative up- and downregulation of cellular pathways. 
Hereby, abundance changes of all factors assigned to 
each pathway are compared to the overall abundance 
distribution and relative values are calculated [21]. Add- 
itionally, significantly up- and downregulated pathways 
were validated by the annotation enrichment tool DAVID, 
which employs a different algorithm for the analysis (data 
not shown) [22]. As expected, pathways altered and 
enriched in the 2-dimensional annotation analysis were 
also identified as enriched by DAVID in all cell lines. 

Our analysis identified a specific pattern of pathways 
that are altered in all compared model trisomic and tetra- 
somic cell lines (Table 1). Among the most upregulated 
categories we identified the endoplasmic reticulum (ER), 
Golgi apparatus, lysosomes and vacuoles, membrane 
metabolism and the MHC protein complex and anti- 
gen processing. In contrast, DNA and RNA metabolic 
pathways - e.g. DNA replication, repair, transcription or 
RNA splicing - were significantly downregulated (Figure 1 A, 
Additional file 1: Figure SIA). Remarkably, although similar 
pathways are up- and downregulated in all aneuploid cell 
lines, we noticed that the specific factors that significantly 
contribute to the differentiated regulation are variable. 
For example, the members of the proton oligopeptide 
cotransporter family SCL15 (Gene Ontology Cellular 
Component "membrane, lysosome, lytic vacuole") show 
variable expression levels in the analyzed aneuploid cell 
lines: SLC15A3 is upregulated in all of the HCT116 de- 
rived aneuploid cell Unes, but not in the RPEl derived 
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Table 1 Pathway classes recurrently and significantly altered in aneuploid cell lines 
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Asterisks indicate the cell lines expressing H2B-GFP. 
/ = annotation not significantly altered. 
+= annotation category upregulated. 
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aneuploid cell lines, whereas SLC15A4 is upregulated 
only in HCT116 H2B-GFP 5/4 and SLC15A1 is upregu- 
lated in three out of seven aneuploid cell lines derived 
from HCT116 and RPEl (Additional file 2: Table SI). 

Interestingly, DLDl with trisomy of chromosome 7 
shows only partial overlap with the other aneuploid cell 
lines: pathways related to splicing and to DNA and RNA 
metabolism were similarly downregulated, and vacuole 
was upregulated, whereas additional changes in the tran- 
scription profile differed from the other trisomic and 
tetrasomic cell lines (Additional file 1: Figure SIB, Table 1). 
We also determined the response to aneuploidy in three 
clones with trisomy of chromosome 8 (here labeled as 
HE35 8/3_l, 2, and 3) derived from human primary 
embryonic fibroblasts HE35 by microcell-mediated chro- 
mosome transfer [17]. Using the 2-D annotation analysis 
algorithm, we observed similar changes in pathway regu- 
lations in the HE35-derived trisomies as in the other 
tri- and tetrasomic cell lines (Figure IB, Additional file 1: 
Figure SIC, Table 1). Taken together, the results indicate 
that all the analyzed human cells with de novo created 
trisomies and tetrasomies elicit a nearly identical pattern 
of pathway changes regardless of the identity of the super- 
numerary chromosome and the cell line. 

Complex aneuploidy triggers the same pathway changes 
as low-complexity aneuploidy 

Next, we asked whether similar pattern of changes in 
gene expression could be observed in cells with complex 
aneuploidy. This type of karyotypic changes might more 
closely resemble the situation observed in tumors that 



often harbor multiple changes in chromosome numbers 
and structures, including deletions, translocations and 
amplifications [23]. To this end, we used cell lines derived 
from HCT116 that underwent a transient tetraploidy 
induced by cytokinesis failure (Figure 2A). As tetraploidy 
leads to catastrophic mitosis and chromosomal instability, 
the majority of cells die soon after the tetraploidization 
[9,24], and the chromosome numbers of the few survivors 
are remarkably altered (Figure 2B). We analyzed the tran- 
scriptional changes identified in two clonal survivors, 
cell lines HPTl and HPT2, by calculating the fold 
change of the transcripts in comparison to the original 
HCT116. This first global analysis of the transcriptional 
response in human cells with model complex aneuploidy 
reproduced the pattern of pathway changes observed in 
the tri- and tetrasomic cell lines (Figure 2C,D, Additional 
file 1: Figure SID, E, Table 1, Additional file 3). Im- 
portantly, the duration of the cell cycle of the complex 
aneuploid cell lines is indistinguishable from diploids 
(Additional file 1: Figure S2). Previously, it has been 
proposed that the transcriptional changes identified in 
aneuploid cells are caused by the slow growth, which is 
typical for model trisomic cell lines [3,14,15]. Thus, 
our finding indicates that the response to aneuploidy 
in mammalian cells is not always associated with a 
slower progression through the cell cycle, and thereby 
the slower proliferation cannot be the only cause of 
the identified transcriptional changes. 

Distance matrix of the Spearman rank correlation of 
all analyzed cell lines confirmed the high similarity 
among the transcriptional changes triggered by aneuploidy 
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Figure 1 Pathway alterations in tri- and tetrasomic cell lines. 

A. 2-dimensional pathway analysis comparison of HCTl 16 with 
tetrasomy of chromosome 5 and DLDl with trisomy of chromosome 
13. Each dot represents an annotation as defined by GO and KEGG 
databases, color-coding marks related annotations as stated in the 
inset legend. The label "Others" indicates annotations that were not 
generalized in larger categories (see Methods). Examples of 
specific pathways categorized by the color- code are annotated. 

B. 2-dimensional pathway analysis comparison of HE35 (clone 2) 
with trisomy of chromosome 8 and RPE* with trisomy of 
chromosome 21. Asterisks indicate the cell lines with H2B-GFP. 



(Figure 2E), We call the identified common pattern of 
transcriptional changes the Aneuploidy Response Pattern 
(ARP). Our results demonstrate that trisomy can serve as 
a valuable model for analysis of complex aneuploidy. 
Moreover, the identification of ARP suggests that chromo- 
somal imbalance itself might be exploited as a novel target 



for cancer therapy regardless the type of aneuploidy and 
the mechanism by which it was generated. 



Comparison of the aneuploidy response pattern with 
transcriptional responses to stress stimuli 

What triggers the specific transcriptional response to 
aneuploidy in human cells? Current analysis in budding 
yeast suggests that mutations interfering with cell pro- 
liferation result in similar gene expression changes as 
aneuploidy [19]. Alternatively, the transcription changes 
may be triggered by cellular stress caused by the presence 
of an extra chromosome. Thus, we asked whether there 
are stress stimuli that trigger reprogramming of gene 
expression similar to the ARP. To this end, we compared 
the transcription profiles of HCTl 16 cells exposed to 
sub-lethal concentrations of hydrogen peroxide, nitric 
oxide, hydroxyurea, actinomycin D and bafilomycin Al, 
or grown under hypoxic conditions as well as in medium 
with either low or high glucose (for further details see 
Methods). Strikingly, the majority of the stressors trig- 
gered pathway changes that are remarkably different 
from the ARP (Figure 3 A, Additional file 1: Figure S3). 
In contrast, we observed a partial similarity between 
the ARP and the transcriptional changes in HCTl 16 
treated with actinomycin D, a polypeptide antibiotic that 
inhibits the activity of RNA polymerases (Figure 3B,D). 
The 2-D annotation analysis revealed that treatment 
with actinomycin results in downregulation of DNA- and 
RNA- metabolism as well as in upregulation of ER, 
membrane metabolism and lysosome. This might indicate 
that the downregulation of DNA- and RNA metabolism 
in aneuploid cell lines is due to transcriptional inhibition 
and its consequences. However, the overlap includes only 
a few genes with a low Spearman correlation coefficient 
and thus the distance of the Spearman rank correlation is 
higher compared to the distances between the aneuploid 
cell lines (Figure 3D). Remarkably, treatment with 
bafilomycin Al that inhibits vacuolar- ATPases showed 
transcriptional changes nearly identical to the ARP 
(Figure 3C,D). Bafilomycin Al impairs vesicle fusion 
[25] and thus inhibits the final steps of autophagy due 
to the failure of autophagosome-lysosome fusion. As 
in aneuploids, the transcriptional changes include down- 
regulation of DNA and RNA metabolism, whereas mem- 
brane associated annotations are upregulated. In contrast, 
treatment with bafilomycin does not lead to upregulation 
of ER, Golgi or lysosomal pathways (Figure 3C). The simi- 
larity with transcriptional effects of autophagy-inhibiting 
drug suggests that one of the main consequences of aneu- 
ploidy in human cells is the overload of the autophagic 
pathway. Taken together, our findings imply that the tran- 
scriptional changes in response to aneuploidy differ from 
most of the common stress responses, but show shared 
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Figure 2 Complex aneuploidy triggers similar pathway changes as trisomy or tetrasomy. A. Schematic depiction of tlie construction 
of aneuploid cell lines with complex karyotype changes. B. Chromosome numbers in clonal cell lines originating from an unstable tetraploid 
intermediate. The numbers above the plot indicate the median chromosome number. C. 2-dimensional pathway analysis comparison of the complex 
aneuploid cell line HPTl'' and HCTl 16 5/4. D. 2-dimensional comparison of complex aneuploid cell line HPT2* and DLDl 13/3. E. Distance matrix of 
the Spearman rank correlation of all analyzed aneuploid cell lines. Large distances present a low correlation and vice versa. The Spearman 
distance metric is defined as distanceSpearman = (1- correlationSpearman)/2 in the distanceMatrix function as part of the R Class Discovery 
package. The function calculates the distance of the correlation between column vectors, each of which represents different microarray experiments. 
Asterisks indicate the cell lines with H2B-GFP. 
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features with response to conditions limiting autophagy 
and transcription. 

Markers of aneuploidy 

Our results imply the existence of factors with recurrent 
expression changes triggered by aneuploidy that might 
serve as general markers of aneuploidy in human cells. To 
identify genes that are consistently up or downregulated 



in aneuploid cells, we merged all available datasets of 
aneuploid cell lines created in our laboratory. This 
yielded 18 genes whose expression is more than 1.4 fold 
increased and 5 genes whose expression is more than 
1.4 fold decreased in all analyzed aneuploids (Figure 4, 
Table 2). The changes in mRNA levels are consistent 
with proteomics measurements in 7 out of 8 proteins for 
which the data are available [15] (Table 2). Importantly, 
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Figure 4 Expression levels of the putative aneuploidy markers. Heatmap of fold change ratios (aneuploid/diploid) of the 18 upregulated and 
5 downregulated genes, hierarchical clustering of Euclidean distance. Asterisks indicate the cell lines with H2B-GFP. 




expression of only one of these 23 genes (FRY) is similarly 
upregulated by the stress stimuli, suggesting that the 
recurrent expression changes are indeed in response to 
aneuploidy. Using quantitative real-time PGR we validated 
the expression levels of four of the identified upregulated 
genes in eight aneuploid cell lines and corresponding 
diploid controls (RAB27B, COL13A1, HOXB5, GLRX). 
This approach confirmed overexpression of all four genes 
in response to aneuploidy. Interestingly, none of the tested 
transcripts was confirmed in one of the complex aneu- 
ploid cell line, HPTl (Figure 5). These candidate markers 
could potentially enable the discrimination of tumors with 
low frequency of aneuploid cells from tumors with high 
levels of aneuploidy and thus higher malignant potential. 
Further experiments will be required to determine 
whether the candidate markers can be used for clinical 
purposes for example by standard immunohistochemistry. 

Discussion 

Significant recent progress in the analysis of aneuploidy 
in eukaryotes has revealed multiple novel features linked 
to aneuploidy such as growth defects, abnormal protein 



homeostasis and increased genomic instability. These fea- 
tures appear to be conserved among eukaryotes, an obser- 
vation which is further supported by the finding that 
the transcriptional response to aneuploidy in different 
species remarkably correlates [19]. We analyzed in detail 
the transcriptional changes in aneuploid human cells by 
comparing the calculated aneuploid-to-diploid ratio. The 
comparison of eleven aneuploid model cell lines derived 
from four different diploid progenitor cell lines by either 
chromosome transfer of seven different chromosomes or 
by clonal propagation from unstable tetraploid progenitors 
revealed a striking similarity among the cell lines. The 
identified signature of pathways that was altered in all 
analyzed cell lines, which we term the aneuploidy response 
pattern - ARP - is characterized by upregulation of the ER 
and Golgi related pathways, lysosome and lytic vacuoles, 
MHC protein complex and antigen processing, whereas 
DNA and RNA metabolism and ribosome-related path- 
ways were always downregulated (Figures 1, 2, Tables 1, 
2). Our results reveal that the transcriptional changes to 
aneuploidy in human cells are not mediated by a canonical 
p53-dependent stress response, because DLDl and its 
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Table 2 List of consistently up- and downregulated transcripts (cut off 1.4 fold change) 



Gene 



Chromosome 



Function 



Protein 
levels 



Link to cancer? 



9 
10 
11 
12 
13 
14 
15 
16 
17 



PLAUR 
RAB27B 
P4HA2 

FRY 
BDKRBl 
H0XB5 

GDF15 

OASL 
SERPINE2 

IFI44 
AMDHDl 
TMEM171 

GLRX 
TMEM169 

DMBTl 
C0L13A1 

SH2D1B 
MIA 



19 
18 
5 
13 
14 
17 

19 
12 
2 
1 

12 
5 
5 
2 

10 
10 
1 



Plasminogen activator, urokinase receptor 
Secretory GTPase 
Procollagen-proline 
Regulator of actin cytoskeleton 
Bradykinine receptor 
Transcription factor 

Growth differentiation factor 
2-5-oligoadenylate synthetase-like 
Serpine protease inhibitor 
Interferon-induced protein 44 
Protein with amidohydrolase domain 
Transmembrane protein 171 

Glutaredoxin 
Transmembrane protein 169 
Membrane glycoprotein 
Collagen type XIII, alphal 
Signal transduction control 
Melanoma inhibitory activity 



4/5 
4/4 
4/6 
0/3 
NA 
NA 

1/2 
NA 
4/4 
NA 
NA 
NA 
4/4 
NA 
NA 
1/1 
NA 
NA 



Ovarian and colorectal cancer 
Breast cancer 
Metastasis 

Breast and lung cancer 

Leukemia, ovarian carcinoma and 
others 

Ovarian and prostate cancer 
Metastasis 



Pancreatic cancer 



Multiple cancers 



Neuronal tumors 



XYLB 

L0XL3 

MYB/NFIB 
fusion 

EEPD 
ARHGEF39 



3 Xylulokinase homologue NA 

2 Lysyl oxidase homologue NA 

6 Myb - NFIB fusion NA 

7 Endonuclease/exonuclease/phosphatase domain NA 

containing 

9 Rho guanine nucleotide exchange factor 39 NA 



The column "Protein levels" shows how often corresponding changes of the protein levels were Identified by Stingele et al. [15] (number of cell lines with altered 
protein expression/number of cell lines where the protein was quantified; NA - proteomlcs data not available). The rows 1-18 show upregulated genes; the rows 
1- 5 (below) show downregulated genes. 



aneuploid derivatives carry mutant p53. A remarkable 
feature of the identified transcriptional changes is that 
although similar pathways are up- and downregulated 
in individual aneuploids, the individual genes whose 
expression is altered are often variable among the dif- 
ferent cell lines. Furthermore, several pathways were 
up- or downregulated in the analyzed cells that were 
not recurrent and thus cannot be regarded as a general 
consequence of aneuploidy per se. These changes are 
probably specific to the type of extra chromosomes as 
well as to the clonal selection. We also show for the 
first time that the pathway changes in human cells with 
complex, hypotetraploid karyotypes correlate with the 
changes identified in simple trisomies. Thus, the analysis 
of model trisomic cell lines might provide important 
insights into the role of the more complex aneuploidy 
that is frequently found in cancer. 



Two of the analyzed cell lines (trisomy 7 in DLDl and 
one of the three clones of trisomy 8 in HE35) showed 
only a partial aneuploidy response pattern (Additional 
file 1: Figure SIB). There are several possible scenarios 
to explain this finding. First, loss of some chromosome 
parts or a mosaic aneuploid population could lead to a 
weakening of the ARP. However, the transcription levels 
of the transcripts coded on the extra chromosomes were 
elevated in these cell lines as expected, which excludes 
this possibility. Alternatively, features specific to the 
supernumerary chromosome can be responsible for the 
difference, such as in case of chromosome 7. This has 
been previously observed in budding yeast, where some 
disomes, e.g. disome of chromosome 1, do not trigger 
the environmental stress response pattern (ERS) other- 
wise observed in other disomic strains [3]. In support 
of this possibility is the observation that chromosome 7 
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A RAB27B B C0L13A1 

30 T 




HCT116 HCT116* RPE1 RPE1* HCT116 HCT116* RPE1 RPE1* 



C H0XB5 D GLRX 

14 T 




HCT116 HCT116* RPE1 RPE1* HCT116 HCT116* RPE1 RPE1* 



Figure 5 Quantitative real-time PCR analysis of putative aneuploidy marlcers. A. Fold change mRNA expression levels of C0L13A1 
quantified by real-time PCR. B. Fold change mRNA expression levels of RAB27B. C. Fold change mRNA expression levels of H0XB5. D. Fold change 
mRNA expression levels of GLRX. Asterisks indicate the cell lines expressing H2B-GFP. 



is one of the few chromosomes which is more often 
gained than lost in chromosomally unstable cancer cells 
[26]. Thus, an extra copy of chromosome 7 likely repre- 
sents a lesser burden for the cells than most of the 
other chromosomes. Finally, the deviation from the ARP 
might arise due to a mutation occurring during the clonal 
selection after the chromosome transfer. This is likely 
the case in HE35 8/3, where two clones correlate stron- 
ger with other aneuploids than the third trisomy 8 
clone (Figure 2E). Despite these exceptions, our results 
show that aneuploidy triggers a uniform transcriptional 
response in human cells that is independent of the 
identity and quantity of the extra chromosomes and 
their combination as well as of the cell line type. 

The ARP is characterized by strong downregulation of 
DNA and RNA metabolism, which correlates with the 
slow growth observed in nearly all aneuploid cells analyzed 



to date. Indeed, many factors required for DNA replication 
are less abundant, such as the heterohexameric replicative 
helicase MCM (Additional file 2: Table SI), and this may 
partially explain the previously observed slow progression 
through the S-phase [15]. However, the complex aneuploid 
cell lines progress through interphase at a rate closely 
matching the rate of the corresponding diploid cell 
lines and yet the levels of factors involved in replication 
are similarly decreased (Additional file 1: Figure S2). 
Therefore, we conclude that the proliferation impairment 
cannot be the only cause of ARP. 

DNA repair is significantly downregulated, which might 
lead to increased accumulation of DNA damage. The 
abundance changes in the DNA replication and repair 
factors might explain how aneuploidy increases genome 
instability, as has been observed in yeast [4] as well as 
in human aneuploids (S.S., V.P., Z.S., unpublished results). 
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Further research should elucidate what are the causes 
of the consistent downregulation of DNA metabolism 
pathways and what consequences it brings upon aneu- 
ploid cells. 

Among the most upregulated pathways we observed 
the Golgi network, ER related pathways, lysosomes and 
lytic vacuoles, membrane metabolism and the MHC 
protein complex and antigen processing. The function 
of MHC protein complexes is to display fragments of 
proteins from within the cell to immune cells [27]. Cur- 
rently we do not understand the reasons for the increased 
levels of the MHC protein complex transcripts in aneu- 
ploid cells. We hypothesize that increased protein expres- 
sion and degradation elevates the peptide presentation 
by MHC complex, thus elevating the immunogenicity of 
aneuploid cells, similarly as it has been recently observed 
in murine aneuploid tumors [28]. In the future it will be 
important to confirm this observation on the protein level 
and to determine whether the MHC proteins are correctly 
localized and functional. 

Elevated expression of ER related genes might suggest 
an expansion of the ER that was shown in budding 
yeast to occur in order to alleviate ER stress [29]. This 
phenomenon is accompanied by elevated lipid biosyn- 
thesis, but unchanged amounts of ER chaperones [29], 
which closely resembles the transcriptional changes 
observed in aneuploid cell lines. This suggests, together 
with the fact that mouse trisomic MEFs are more sensitive 
to the Hsp90 inhibitor 17-AAG [30], that aneuploid cells 
suffer from protein folding defects, most likely due to the 
saturation of the cellular folding capacity [31]. ER stress 
is often observed in various cancers, where it is usually 
attributed to tumor microenvironments characterized by 
hypoxia, nutrient limitation and low pH [32]. It will be 
interesting to test whether chromosomal copy number 
changes also contribute to the ER stress in cancer cells. 

Previously, it was proposed that the transcriptional 
response to aneuploidy is indicative of stress and slow 
growth [1]. Therefore, we were interested whether 
treatment with stress inducing agents will elicit similar 
response as the presence of an extra chromosome. Since 
chromosome missegregation was shown to trigger oxi- 
dative stress [24,33], we analysed the transcriptional 
response to oxidative stress inducers nitric oxide and 
hydrogen dioxide. Aneuploidy was also recently linked 
to DNA damage and replication stress [34], therefore 
we analyzed the response to the replication inhibitor 
hydroxyurea and transcription inhibition by actinomycin 
D. Further, aneuploidy was shown to trigger proteotoxic 
stress and to activate autophagy [3,15,30,31]. Therefore, 
we compared the response to aneuploidy to the tran- 
scriptional changes elicited by cells grown under hypoxic 
conditions, which is known to result in energy and ER 
stress with subsequent proteotoxic stress. Similarly, we 



tested the transcriptional changes in cells where au- 
tophagy was inhibited, a condition known to trigger 
proteotoxic stress. The comparison with transcription 
profiles of cells grown on high or low glucose was of 
particular interest, because aneuploid cells show changes 
in the metabolic pathways and higher energy demands 
[3,14,15]. Remarkably, we found that none of the stress 
conditions triggers response similar to ARP. 

The only significant exception was the striking similarity 
between the ARP and the transcriptional changes observed 
in cells treated with sub-lethal concentrations of bafilomy- 
cin Al (Figure 3). Treatment with bafilomycin Al leads to 
the accumulation of vesicles and membranes and reduces 
cell proliferation. The similarity of cellular response to 
inhibition of autophagy and to aneuploidy might suggest 
that aneuploidy inhibits autophagy. This however is not 
true as we previously documented by an observation that 
autophagic flux (the dynamic flow through autophagy) is 
in aneuploids as efficient as in diploids [15]. Inhibition 
of basal autophagy by bafilomycin Al treatment causes 
proteotoxic stress [35]. Current data suggest that aneu- 
ploidy leads to proteotoxic stress as well [31]. We propose 
that the similarity between the transcriptional response to 
treatment with bafilomycin Al and the ARP is because 
both reflect transcriptional changes in cells suffering from 
proteotoxic stress. In general, the effect of bafilomycin Al 
further emphasizes the role of autophagy in mammalian 
aneuploids. Aneuploidy leads to activation of autophagy 
in human and mouse cells [15,30] and aneuploid cells 
are more sensitive to the autophagy inhibiting drug 
chloroquine [30]. The consistent upregulation of lysosomes 
and lytic vacuoles appears to be specific for mammalian 
cells as it was not observed in aneuploid yeasts and plants 
[19]. In this context it is interesting that there is nearly no 
correlation of ARP and transcriptional changes upon low 
glucose condition, which is known to trigger autophagy 
(Additional file 1: Figure S3). We believe that this is 
because autophagy in aneuploids is not activated by energy 
or nutrition deprivation, but by different, so far unidenti- 
fied triggers [36]. Partial overlap of ARP was also observed 
with the response to treatment with the drug actinomycin 
D that represses transcription by RNA polymerase and, 
among other effects, leads to an imbalance of ribosomal 
subunits [37]. Transcription-related pathways are also 
downregulated in aneuploids, suggesting the possibility 
that some of the phenotypic features of aneuploid cells 
might be caused by transcriptional deficiency. 

Our work has also allowed the identification of several 
genes that were consistently up- or downregulated in an- 
euploid cell lines, but not in the HCT116 cells under 
stress stimuli. Several of the identified factors have been 
previously linked to cancer (Table 2). For example the 
plasminogen activator and urokinase receptor PLAUR/ 
UPAR is frequently overexpressed in ovarian and 
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colorectal cancers, where it facilitates cell motility and 
metastatic potential and emerges as a marker of malig- 
nancy [38]. P4HA2 is a precursor of collagen and a metas- 
tasis marker [39] and HOXB5 is a homeobox 
transcription factor whose overexpression has been 
found in several malignant tumors [40]. Transcription of 
the TGFp superfamily cytokine GDF15 (growth differenti- 
ation factor 15) is increased in multiple solid tumors, 
where it promotes tumor growth and confers resist- 
ance to several drugs such as bortezomib. Recent studies 
proposed GDF15 as a prognostic marker for ovarian and 
pancreatic carcinoma [41], two types of malignancies 
with frequent aneuploidy and chromosomal instability 
[9]. Glutaredoxin (GLRX) protects cells from oxidative 
stress and serves as a potential marker of malignancy 
in pancreatic carcinoma. Interestingly, some of the fac- 
tors such as a xylulokinase homologue (XYLB), the 2'-5'- 
oligoadenylate synthetase like (OASL) and the interferon- 
induced protein (IFI44) have not been so far linked to can- 
cer and may represent markers of aneuploidy that are not 
linked to malignancy. 

Aneuploidy is associated with poor prognosis and in- 
creased drug resistance in tumors [42]. This suggests an 
interesting possibility that some of the genes that are 
expressed at higher levels in response to aneuploidy 
in non-transformed cells may subsequently enhance 
the malignant potential of these cells, whereas higher 
expression of others might represent a barrier for car- 
cinogenesis. The expression of the marker genes, if 
confirmed in a wide range of aneuploid cell lines and in 
cancer cells, might provide an excellent tool for cancer 
biology and treatment, because it may allow distinguish 
tumors with high aneuploidy from predominantly diploid 
tumors. Moreover, pathways that are activated or inhibited 
by aneuploidy may serve as novel targets for cancer 
treatment. This has been recently demonstrated by the 
fact that aneuploid cells are more sensitive to drugs 
inhibiting HSP90, autophagy inhibitors and inhibitors 
of the AMPK kinase, which strongly correlates with the 
identified increased requirements for autophagy, protein 
folding and energy metabolism in aneuploid mammalian 
cells ([15,30] and this work). Further research should 
elucidate the efficacy of this approach. Our work for the 
first time identifies global changes in a broad spectrum of 
human aneuploid cell lines and may therefore help to 
generate new hypothesis for cancer treatments. 

Conclusions 

In this study we identified a transcriptional aneuploid re- 
sponse pattern (ARP), a set of transcriptional changes, 
common in a broad range of human aneuploid model 
cell lines. This, for the first time, shows that complex 
aneuploidy, which is frequently found in cancer cells, 
exhibits the same transcriptional pathway changes as 



simple trisomy and tetrasomy. A general response to 
aneuploidy, as identified here in a variety of aneuploid 
model cell lines, might serve as a novel therapeutic 
target in cancer treatment. Further, we found 23 genes 
consistently deregulated in our model aneuploid cell 
lines. A confirmation of these markers in aneuploid 
cells might open new strategies for identifying aneu- 
ploid tumors. Since the fraction of aneuploid cells in 
tumor correlates with malignant potential and poor 
prognosis in cancer, simple and reliable biomarkers for 
aneuploidy may help for determining the appropriate 
cancer therapy. 

Methods 

Cell lines 

Following model aneuploid cell lines were used for data 
analysis: parental cell line HCT116 (human colon car- 
cinoma cell line): HCT116 3/3 (trisomy 3), HCT116 5/4 
(tetrasomy 5) [43]; parental cell line HCT116 H2B-GFP: 
HCT116 5/4 (tetrasomy 5) [15], HPTl, HPT2 (hyperte- 
traploids with complex karyotypes, A.Y.K. unpublished 
data); parental cell line RPEl (human retinal pigment 
epithelial cell line, hTERT immortalized): RPE15/3 12/3 
(trisomy 5, 12); parental cell line RPEl H2B-GFP: RPEl 
21/3 (trisomy 21) [15]; parental cell line DLDl (human 
colon adenocarcinoma cell line): DLDl 3/3 (trisomy 3), 
DLDl 7/3 (trisomy 7), DLDl 13/3 (trisomy 13) [13], par- 
ental cell line HE35 (human embryonic cell line): HE35 
8/3, clones 1-3 [17]. 

Cell cultures 

The HCT116- and RPEl- derived tri- and tetrasomic cell 
lines have been described previously [15]. The post- tetra- 
ploid cell lines HPTl and HPT2 were generated by expan- 
sion of individual tetraploids formed by induced inhibition 
of cytokinesis through dihydrocytochalasin treatment. The 
DNA content was determined by flow cytometry, standard 
karyotyping, chromosome painting and array comparative 
genomic hybridization as in [15]. Cells were maintained 
in 10 cm dishes in Dulbeccos Modified Eagle Medium 
(DMEM) at 37°C and 5% C02. Culture medium was 
supplemented with growth factors from 5% Fetal Calf 
Serum (FCS) and 1% Penicillin/ Streptavidin to avoid 
bacterial growth and contamination. Specific antibiotics 
were added to the medium when necessary to maintain 
the supernumerary chromosomes. Cell lines were cultured 
at 70- 90% confluence for a maximum of 10 passages. 

Material for transcription analysis was obtained at 
the earliest possible time point, after sufficient amount 
of cells was achieved, approximately 25 generations 
after the chromosome transfer or tetraploidization, re- 
spectively. Similar approach was taken for the analysis 
of trisomy 8 HE35-derived clones [17]. The mRNA of 
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DLDl -derived aneuploids was isolated from multiple 
different passages [13]. 

Cell growth analysis 

Freshly cultured cells that carry H2B-GFP to fluorescently 
label the chromosomes were seeded 24 h before the ex- 
periment. Time laps movies were taken by imaging 
asynchronous cells in a 10 min or 4 min interval for 
72 or 48 h, respectively. The time in interphase was 
measured as the time from nuclear envelope reformation 
to nuclear envelope breal< down. 

Quantitative real-time PGR 

Total RNA was extracted with the RNeasy Mini Kit 
(Qiagen), treated with DNAse (recombinant DNase, Roche) 
and subsequently transcribed into cDNA (Transcriptor 
First Strand cDNA Synthesis Kit, Roche Diagnostics). 
Specific primers were designed using PrimerBlast; the 
sequences are listed below. Quantitative PGR was 
conducted using the Light Cycler 480 System (Roche 
Diagnostics) with the KAPA SYBR FAST master mix 
optimized for Roche Light Cycler 480 (KapaBiosystems). 
Absolute quantification with an external standard was 
performed and negative non-template controls were 
tested in all experiments. The specificity of the primer 
product amplification was confirmed in each run by 
melting curve analysis. mRNA expression was normalized 
to the control gene coding for ribosomal protein L30 
(RPL30) and fold change to corresponding diploid mRNA 
expression was calculated. Primers: COL13A1 - forward: 
GGGGGAAGCAGGACTAGATG, reverse: CCTGAAGCT 
CCGGGTAGTQ RAB27B - forward: TGCGGGACAAGA 
GCGGTTCCG, reverse: GCCAGTTCCCGAGCTTGCCG 
TT, HOXB5 - forward: TCCACAAATCAAGCCCTCCA, 
reverse: GTCCGGGCCATTTGGATAAQ GLRX - forward: 
AACGGTGCCTCGAGTCTTTA, reverse: CCTATGAGATC 
TGTGGTTACTGG 

mRNA expression analysis by microarrays 

Genome-wide expression profiling of HCT116- and RPEl- 
derived aneuploid cells lines was conducted in three 
replicates by IMGM laboratories GmbH (Martinsried, 
Germany) as previously described [15]. cDNA was hybrid- 
ized on Agilent Whole Human Genome Oligo microarrays 
(4x44K format) for HCT116 diploid and HCT116 5/4, 
or Agilent SurePrint G3 Human GE microarrays (8x60K) 
for the other cell lines according to a One-Color based 
hybridization protocol. Raw data was background normal- 
ized. The data has been deposited in the NCBI Gene 
Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/ 
geo/) accessible through GEO Series accession number 
GSE47830 and GSE47836. 

Microarray data of trisomic and diploid colorectal can- 
cer cell lines DLDl [13] were kindly provided by Thomas 



Ried (National Institutes of Health, Bethesda, Maryland, 
USA). All other mRNA expression data were obtained 
from NCBIs Gene Expression Omnibus (http://www.ncbi. 
nlm.nih.gov/geo/). HCT116 cell line grown under high 
and low glucose conditions: accession number GSE31084; 
HCT116 after stress treatment: accession number GSE3176; 
HCT116 cell line treated with actinomycin D: with 
GSE12459, human embryonic diploid and trisomic 
cells: accession number GSE28076. 

Microarray normalization 

Bioinformatics analysis of the microarray data was per- 
formed using Perseus (1.2.6.16) as part of the MaxQuant 
Software Package [21] and R on the free open source 
integrated development environment R Studio. The gBG 
SubSignal (green Background Subtracted Signal) from 
background-normalized cDNA microarray data was used 
for further normalization. The background-subtracted raw 
intensities were log2 transformed and global normal- 
ization of the log transformed raw data was performed by 
subtracting the median of the overall signal intensities for 
one experiment from each signal in this experiment. Probe 
sets for one gene were summarized by taking the median. 
The median of replicative probes signal intensities was 
calculated for each cell line. For comparison of each 
aneuploid cell Une with the corresponding diploid cell 
line gene, expression fold change ratios were calculated. 

Data analysis 

The Students ^-test was performed to verify the statis- 
tical significance of the fold change in mRNA expression 
between the signal intensities of aneuploid and diploid 
cell line. To correct the test statistics for multiple com- 
parisons, false discovery rate control of the p-values was 
applied. Both local and frequent FDR were calculated 
with the "fdrtool" package [44]. For further analysis, a 
fold change cut off 1.4 in mRNA expression was applied. 
To investigate the correlation between the cell lines, a 
distance matrix of the Spearman rank correlation coeffi- 
cient was calculated in R using the "ClassDiscovery" 
package in the OOMPA library and "mclust" package. 
Distance of the correlation was visualized in a colored 
matrix in the "spatstat" package. To investigate a global 
response to aneuploidy on gene level, data sets were 
filtered for those genes more than 1.4 fold up- or down- 
regulated in all HCT116 and RPEl derived cell lines. 
The resulting gene lists were analyzed for their function 
and relation using the DAVID functional annotation tool 
[22], KEGG pathway and the NCBI Gene database. 

Pathway enrichment analysis 

Pathway enrichment analysis was conducted in Perseus 
using the 2-dimensional annotation enrichment tool [21]. 
Thereby, enriched GO and KEGG annotations were 
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identified by testing whether genes in an annotated cat- 
egory have a significant preference to be altered com- 
pared to the global fold change data distribution. 
Significance was tested by a ^-test followed by a false 
discovery rate correction in the Benjamini-Hochberg 
procedure (FDR cut off 0.02). Expression values for sig- 
nificant enriched annotations were summarized in an 
annotation score from -1 to 1 to represent levels of up- 
or downregulation. Significant and enriched annotations 
were summarized in larger, general categories for 
visualization purposes. 2-dimensional annotations en- 
richments of two cell lines were plotted in R. Highly up- 
or downregulated pathways were cross validated by sub- 
mitting the fold change cut off datasets to the annotation 
enrichment tool DAVID [22]. 

Statistical analysis 

Commonly used R packages were "lattice" [45], "genefilter" 
[46], "fdrtool" [44] and "calibrate" [47]. In addition, Perseus 
(1.2.6.16) as part of the MaxQuant Software Package [21] 
was used for microarray analysis. All statistically evaluated 
experiments were performed in at least three independent 
biological replicates. The final plots were prepared in 
GraphPad Prism 5 software or R. 
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